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ABSTRACT 


The Space Kinetic Impact and Debris Branch 
has begun an ambitious program to con- 
struct a fully analytical model of the 
breakup of a satellite under hypervelocity 
impact. In order to provide empirical 
data with which to substantiate the model, 
debris from hypervelocity experiments 
conducted in a controlled laboratory envi- 
ronment has been characterized to provide 
information on its mass, velocity, and 
ballistic coefficient distributions. Data 
on the debris has been collected in one 
master data file, and a simple FORTRAN 
program allows users to describe the de- 
bris from any subset of these experiments 
that may be of interest to them. A sta- 
tistical analysis has been performed, 
allowing users to determine the precision 
of the velocity measurements for the data. 
Attempts are being made to include and 
correlate other laboratory data, as well 
as those data obtained from the explosion 
or collision of spacecraft in low earth 
orbit . 


INTRODUCTION 


Characterization of debris from hypervel- 
ocity impact events is an important pre- 
requisite for analytical modelling of 
those events. One feature of a useful 
model would be its ability to predict the 
characteristics of the debris cloud pro- 
duced by the impact, and this feature 
cannot be evaluated without a body of 
empirical data with which to compare the 
analytical predictions. The Space Kinetic 
Impact and Debris Branch (SKID) has con- 
ducted a debris characterization program, 
using hypervelocity impact debris produced 
in antisatellite experiments. Our ap- 
proach includes some novel ideas from 
which the orbital debris community may 
benefit . 


WEIGHING AND CATEGORIZING OF DEBRIS 


SKID has characterized debris from hyperv- 
elocity shots conducted at the Naval Re- 
search Laboratory (NRL) and at the Uni- 
versity of Dayton Research Institute (UDR- 
I) . Most debris comes in plastic bags of 
two sizes, large and small. The large 
bags contain debris that has been swept 
from the floor of the impact chamber, and 
include a considerable amount of extrane- 
ous material. The smaller bags contain 
debris recovered from the inside of the 
impact target, and do not contain as much 
unwanted material. 

The characterization begins with the sort- 
ing out of the metallic debris, from the 
non-metallic debris; this is done by hand. 
This provides another opportunity to re- 
move extraneous material. The metallic 
debris is then weighed on the microbalance 
at Phillips Laboratory's metrology facili- 
ty. This microbalance has an accuracy of 
0.004 g. Debris particles with mass less 
than 10 times the accuracy of the scale 
are not weighed. The mass of the parti- 
cle, as well as all pertinent information 
about the shot from which it came (shot 
number, impact angle, nominal projectile 
velocity, projectile material), is record- 
ed in a master data file. 

Some debris was collected in catcher mate- 
rial during the test, for the purpose of 
measuring its velocity. This debris is 
removed from the catcher material with 
tweezers, and the distance that it trav- 
elled into the catcher material is mea- 
sured. This distance is used to determine 
the velocity at which the particle had 
travelled, using methods to be discussed 
later in this paper. The mass, velocity, 
and shot information are then included in 
the master data file. 
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CALCULATION OF DISTRIBUTION CURVES 


Calculation of mass distribution curves is 
performed by the FORTRAN-77 code DEBRIS. F. 
The program selects the data of interest 
to the user, and then generates a mass 
distribution plot, which is written to a 
file called DISTRIB.DAT. The data can be 
plotted in two ways. In the first ap- 
proach, the dependent variable is mass, 
and the independent variable is the total 
number of debris particles larger than 
that mass. Such a plotting method is 
often used in the space debris community, 
for example, in References 1 and 2. In 
the second approach, the dependent vari- 
able is still mass, but the particles are 
separated into "bins.” For example, a bin 
might contain all particles with mass 
between 0.1 g and 0.3 g. If 12 particles 
are found to lie in this range, the re- 
sulting data point would be (0.2, 12). 

This method of plotting is similar to the 
approach used by the debris hazard model- 
ling computer program IMPACT [3] . 

DEBRIS. F also fits a straight line and a 
parabola to the curve, and performs an F- 
test (as described in Reference 4) to 
determine if the higher-order fit is need- 
ed. The least-squares coefficients and 
the F-test results are written to a file 
called FITS . DAT. The source code can also 
generate velocity and ballistic coeffi- 
cient distributions. DEBRIS. F, and the 
master data file DEBDAT, are installed on 
the CRAY-2 supercomputer at Kirtland AFB , 
whence they can be easily downloaded onto 
floppy disks. 

It has been noted by McKnight c5] that a 
straight-line fit to a mass distribution 
does not make physical sense: if the line 

is extrapolated to the left, it is seen to 
imply the presence of an infinite number 
of infinitesimal particles. It makes more 
sense for the curve to have a "knee." 
Such a "knee" can be introduced with a 
bilinear least-squares fit, which has been 
derived as a part of this effort, and is 
available alonq with the linear and para- 
bolic fits as part of the DEBRIS. F code. 

Consider the situation shown in Figure 1: 
a bent straight line is to be fit to a set 
of data. The line is defined by four 
parameters: the intercept p 0 , the two 

slopes P y and P 2 , and the x-coordinate of 
the break point, x 0 . Consider the case in 
which x 0 is given. A least-squares fit is 
obtained by minimizing the residual sum of 
squares SS ( res), given by 



SS(res) = £ [yi-(P 0 + Pi x i)] 2 

i-1 

+ [ (yi“Po - Pl*o) - p2 


( 1 ) 


where x < x 0 for data points 1 through n s , 
and n is the total number of data points. 
This T is minimized by setting the partial 
derivatives of SS(res) with respect to the 
parameters equal to zero: 


a SS(res) = a 5 5 (res ) 

ap 0 “ a Pi 

a 55 (res) = , 

ap 2 


( 2 ) 


Substituting equation 1 into equations 2 
leads to a system of three linear equa- 
tions in three unknowns: 


[a] 


Po 

Pi 

p 2 


[y] 


(3) 


where [cr] is a symmetric 3x3 matrix: 


®11 Tlrr 


( 4 ) 


n, 

a i2 = n N X o + 5^ x i 
i-l 


(5) 
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n N X 0 


«13 = -n„x o + £ x i 

i-n,-*-! 


«22 = >V*0 2 + £ X? 


Figure 2 shows a mass distribution from 
(6) DEBRIS. F, with the linear and parabolic 

fits employed. Data from all shots that 
were characterized was combined into one 
file, which is plotted here. The x-coor- 
dinate is a chosen mass M. The y-coordi- 
d) nate is the number of particles observed 

to have a mass greater than M. The F-test 
performed by the program indicates that 


T3 r 

«23 = ~n N 4 + x 0 £ x 1 

i m n B *l 


= £ Uj-Xo) 1 


i m n M * 1 


the vector [y] is defined as 


[y] 




£ <x' J -x 0 )y i 


i-n 0 + 1 



Figure 2. Mass distribution curve. 


( 10 ) 


the accuracy gained by using the parabolic 
fit, as opposed to the linear fit, is 
statistically significant in this case. 


and n N is the number of data points for 
which x > x Q . This system can now be 
solved for P Qf 0 1# and /? 2 . Note that this 
solution requires a "guess” for Xg, but not 
for the y-value of the break point. The 
FORTRAN code DEBRIS. F performs 1001 bent- 
line fits to a given data set, using 1001 
evenly-spaced guesses for x 0 . The fit with 
the best correlation coefficient is re- 
tained. 

This bilinear fit is particularly useful 
for applying the bilinear exponential 
model of Grady and Kipp c63 . Reference 6 
proposes an equation of the form 

y = + A 2 e‘^ x dU 


This equation cannot be fit to a set of 
data using simple least-squares means; a 
more sophisticated, iterative approach is 
required. However, equation 11 approxi- 
mates the bent-line fit if the x-axis is a 
linear scale, and the y-axis is a loga- 
rithmic scale. Thus, the bent-line fit 
can be applied to the model of Reference 
6 . 


Figure 3 shows the same data, separated 
into bins of uniform width on a logarith- 
mic scale. (Bins can also be of uniform 
width on a linear scale.) The x-coordi- 
nate of each data point corresponds to the 
center of the bin. Linear and parabolic 
fits are also available with this type of 
plot . 


# pirtlclu In bln 



Bin Mean, grams 


— Raw Data ~ Llnaar Fit -*- Parabolic Fit 


Figure 3. Sample plot with data separated 
into bins 


Figure 4 shows the same data, plotted with 
the bilinear fit routine. At first 
glance, the fit does not appear to be very 
good. However, a closer look shows that 
the linear scale on the x-axis leads to a 


744 




concentration of data points near the 
left-hand side of the graph. Once these 
points are recognized, the fit is seen to 
be quite good. 

Figure 4 also shows that, when the distri- 
bution data is plotted with a linear x- 
axis, the resulting curve is almost verti- 
cal at low values of x. For some test 
data, this steepness caused numerical 
stability problems in the bilinear fit 
routine, which resulted in negative values 
of R 2 for some or all of the trial fits. 
These stability problems were decreased, 
but not eliminated, by restricting the 
range of the trial values of x Q . The 
problem was solved by implementing a com- 
pletely pivoted Gaussian elimination rou- 
tine to solve equation 3. DEBRIS. F had 



Figure 4. Mass distribution plot 
with bilinear fit. 


Figures 5 and 6 show individual plots for 
all tests, using the plotting formats of 
Figures 2 and 3. Figure 6 used a linear 
vertical axis, due to the presence of some 
zero values in the data. This format also 
emphasizes the differences in the data 
from the various tests. The most inter- 
esting result is the considerable varia- 
tion among the curves shown in these two 
figures. These shots are actually quite 
similar; yet, their debris distributions 
are quite different. It should be noted 
that, while no debris particles of mass 
less than 0.04 g were weighed, this limi- 
tation would not change the values given 
in Figures 5 and 6; it would merely allow 
the curves to be extended to the left. 
Thus, the variability shown in these fig- 
ures is real . 



Figure 5. Mass distribution plots 
for individual tests. 
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8 hot 0 2-691 Shot 0 2-692 Shot 0 2-696 Shot 0 2-694 

Shot 0 2-696 Shot 0 2-6*6 Shot 0 2-697 Shot 0 2-896 

Figure 6. Individual mass distribution 
plots, using bin format. 

The flexibility of this approach would 
also be useful in studies such as Refer- 
ence 2. McKnight and Brechin determined 
the linear fit coefficients for the mass 
distributions of several impact events, 
and fit a curve to these coefficients in 
an effort to create a master equation to 
describe the mass distributions of a wide 
range of impact events. Using the SKID 
method of cataloguing the debris, a large 
body of data could be built up and easily 
analyzed. Effects of varying parameters 
(such as projectile mass, projectile ve- 
locity, projectile material, etcetera) on 
the linear or parabolic fit parameters 
could be observed. 


STATISTICAL ANALYSIS OF VELOCITY 
MEASUREMENTS 


Velocity measurements are taken using the 
results of Malick's empirical study cn . 
Malick fired test projectiles of four 
different materials at wallboard, using a 


745 




range of known velocities, masses, impact 
angles (see Figure 7) , and presented par- 
ticle areas. 

Malick then fit a curve to his data. The 



Figure 7. Definition of impact angle 0 
( AABC) ’ 


resulting equation was: 


1537 f y t 0 - 8091 
Bl°‘ 3336 (COS 0) 0 - 5419 


( 12 ) 



Figure 8 . Power-law fit to Malick's data 




41 . 2 


(13) 


where p is the projectile density in lb/- 
ft 3 . This curve fits the common logarithms 
of Malick* s data with a correlation coef- 
ficient of R 2 =0 . 9996 . Figure 8 shows a 95% 
prediction confidence interval drawn about 
this curve, illustrating the reliability 
of interpolations of f v made on the basis 
of it. Incorporating equation (13) into 
equation (12) gives 


where V is the particle velocity in ft/- 
sec, f is a function of the projectile 
density, t is the thickness of particle 
board penetrated, in inches, m is the 
weight of the particle, in grains, and 0 
is the impact angle. The multiple corre- 
lation coefficient for the fit of this 
curve to the data is 0.88. Equation (12) 
applies only to the case of approximately 
cylindrical projectiles, for which 

A » c/r? 2/3 (13) 


6 . 33 2 xlO 4 t°- 8091 

/7? 0.3336 p 0.6 (cos 0)0.5419 


It is important to know the reliability of 
velocity measurements made with Malick' s 
equation. A prediction confidence inter- 
val for these measurements can be con- 
structed, using the approach found in 
Bancroft and Han. c83 

A multiple linear regression used to ob- 
tain equation (14) would begin with an 
equation of the form 


where A is the average presented area of 
the particle, and c is a constant. 

Figure 8 shows the results of fitting a 
power-law curve to Malick' s data for f y as 
a function of density. The resulting 
curve is: 


log V f i = b 0 + b 1 logp 1 

+ b 2 log t d + b 3 logm i (15) 

+ Jb 4 log (cosBj) 

where V*. is the least-squares approxima- 
tion to the actual velocity V. . For conve- 
nience, rewrite equation (15) as 


log V / i = b 0 + b lXli + b 2 x 2i 
+ i> 3 x 3J + b t x ti 


(16) 
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where the following abbreviations are 
used : 


logp = x 1 

( 17 ) 

log t = x 2 

( 18 ) 

log/n = x 3 

( 19 ) 

log (cos0) = x 4 

( 20 ) 


V'i ± t (</2 ,, p) s/l«ilxl ( 25 ) 


Thus, it is necessary to calculate the 
estimator S 2 and the matrix [C] from Mal- 
ick's original data, after which the de- 
sired confidence intervals can be calcu- 
lated by plugging the vector {X} into the 
appropriate equation. 

Converting equation (14) to the MKS system 
gives 

where V is now measured in m/s ec, p is in 
kg/m 3 , t is in meters, and m is in kilo- 


Now, define the matrix [C] to be the in- 
verse of the matrix [D], for which 


D kl “ E 

i- 1 


where all summations are taken over all n 
data points { V, , x,. , x,. ,x 3i , x 4i } , and x o is 
defined as 1. Furthermore, define the 
quantity S 2 : 


1.108X10 3 t 0 - 8091 (26) 

n? 0.3336 p 0.6 (cos 0)0.5419 

grams. With all data available from Ref- 
erence 7 converted to MKS, the matrix [C] 
was calculated, and, using equation (26) 
as the fit to the data, the value of S was 
determined to be 

5=0. 144 < 27 > 


5 2 = Y (logVj -log V* ± ) 2 ( 22 ) 

n~P H 


where p is the number of degrees of free- 
dom in the model. For this case, p=5; 
there are four coefficients in the re- 
gression equation, plus one constant. S 
is an unbiased estimator of the variance 
of the errors of the data relative to the 
least-squares fit, assuming those errors 
to be normally distributed. Under this 
assumption, a 100(l-a)% confidence inter- 
val for V. is given by 

v'l ± t,./ 2 . n - P) Sv/W T [c]bri 


where 

{*} = ( 1 x x x 2 x 3 x 4 ) T 


( 23 ) 


( 24 ) 


and t z is the 100(l-a/2)% point of 
Student's t-distribution, with n-p degrees 
of freedom. Tables of this distribution 
are readily available in most probability 
texts. A 100(l-a)% prediction interval is 
given by 


Two facts about the data of Reference 7 
are worth noting. First, the majority of 
shots were fired at impact angles of zero 
degrees, leaving relatively little infor- 
mation on the variation of V with 0 . This 
is reflected in the large value of the 
coefficient C 55 in the matrix [C], which 
indicates that the width of the confidence 
interval is very sensitive to 0 . Second, 
the highest impact velocity recorded for 
valid data is about 9000 ft/sec, or about 
2.7 km/sec. Malick performed some shots 
at higher speeds, but considered the re- 
sulting data invalid because of excessive 
breakup or deformation of the projectiles. 
Thus, the average velocity for Malick' s 
shots is probably 1.0 - 1.5 km/sec. Fig- 
ure 9, a schematic of the confidence and 
prediction bands about a simple least- 
squares straight-line fit, shows why this 
is important. The confidence band, which 
is described by equation 23, is tightest 
at the mean x-value of the data. Its 
width varies as the square of the distance 
from this point along the x-axis. Because 
of the linear nature of the fit, one can 
make a similar statement regarding the 
variation with distance from the mean y- 
value of the data. For the same confi- 
dence level, the prediction band, des- 
cribed by equation 25, behaves in a simi- 
lar manner, but is much wider. The geo- 
metrical interpretation is more complex 
for a multiple regression. However, 
one can say that Malick 1 s fit is probably 
most accurate for velocities of approxi- 
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mately 1.0 - 1.5 km/sec, and that the ac- 
curacy of this method decays for values 
significantly outside of this range. 
Thus, it is important to have a quanti- 



prediction intervals. 


CONCLUSION 

Since the appropriate statistical parame- 
ters have been calculated, one can now use 
equations (23) and (25) to determine the 
accuracy of any velocity prediction made 
using Malick's method. 

The method of creating debris character- 
ization plots is very versatile. The data 
file collected for this program can be 
added to at any time, and the computer 
code that reads it is simple and easily 
modified. This flexible approach will 
allow this database to be expanded and 
restudied as new issues arise in debris 
characterization . 
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